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Оцінювання прогнозів зсувних процесів 
Південного берега Криму за допомогою 
аналізу часових рядів (АРПКС) 


Досліджено загальні закономірності побудови авторегресійних моделей зсувного процесу на прикладі 
Південного берега Криму. Проведено ідентифікацію АРПКС-моделей та розрахунки оцінок параметрів 
цих моделей. Обрано критерії оцінювання моделей та оцінювання прогнозів, за якими виділяється модель 
серед інших, та побудовано прогноз на 10 років. 


Вступ 


Задача аналізу, моделювання та прогнозування зсувних процесів Південного бе- 
рега Криму обумовлена присутністю в даній системі чинника ризику непередбачених 
ситуацій та ускладнена невизначеністю зміни кліматичних, екологічних, геологічних 
та техногенних процесів, які негативно впливають на господарські об'єкти та інженер- 
но-будівельні конструкції, на погіршення їх надійності 1 довговічності, а також при- 
зводять до їх руйнації (| 1. 

Таким чином, стає актуальною проблема системного аналізу динаміки зсувних 
процесів, а також оцінювання прогнозів щодо їх активізації, можливих збитків та обся- 
гу коштів, які знадобляться для покриття матеріальних збитків від руйнації об'єктів 
господарювання внаслідок зсувів. 

Отже, стає необхідною розробка нових моделей, методик, алгоритмів, інтелек- 
туальних систем, які б допомагали особі, що приймає рішення, у виробленні 1 прийнятті 
управлінських рішень в умовах невизначеності. 

У наукових дослідженнях, що присвячені прогнозуванню зсувних процесів Пів- 
денного берега Криму, слабкою ланкою є спрямованість лише на спостереження, мо- 
ніторинг, картування та експертну оцінку фахівця, відсутність системного підходу, об- 
межене використання сучасних інформаційних технологій та інтелектуальних систем 
прийняття рішень. 

Метою даної роботи є розвиток методики прогнозування зсувних процесів Пів- 
денного берега Криму з використанням регресійного аналізу, задля чого виконується 
дослідження і аналіз факторів. Основним критерієм при побудові моделей обрано адек- 
ватність спостереженням, точність прогнозу та стійкість моделі. 

Визначена мета статті зумовила необхідність вирішення таких задач: 

- аналіз факторів та змінних, які описують зсувні процеси; 

- розробка математичних моделей - авторегресії, АРПКС, їх аналіз та оцінку; 

- побудова прогнозу зсувних процесів Південного берега Криму. 
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У рамках даної роботи розглядається задача оцінювання оперативного прогнозу 
зсувних процесів Південного берега Криму, що є складовою задачі управління екзоген- 
ними процесами 1 пов'язана з мінімізацією вимушених витрат на відновлення ушко- 
джених об'єктів. Вона має складний характер внаслідок недостатньої формалізації та 
структурованості, залежності від поточної ситуації, високого рівня невизначеності, об- 
меженості у часі для прийняття рішення. 


Методика Бокса-Дженкінса для побудови 
математичної моделі 


Методика побудови математичної моделі процесу у вигляді різницевого рівнян- 
ня складається з наступних кроків |21: 

- обчислення та аналіз автокореляційної та часткової автокореляційної функцій 
(АКФ та ЧАКФ); 

- оцінювання коефіцієнтів математичної моделі вибраної структури (чи багатьох 
структур) за допомогою вибраного методу оцінювання (найчастіше це МНК або його 
модифікації); 

- діагностика отриманих моделей та вибір з них моделі, адекватної процесу, 
що моделюється. 

Розглянемо детальніше згадані етапи побудови моделі. 

На першому етапі виконують нормування та візуальну перевірку експеримен- 
тальних (статистичних) даних і за необхідності коригують їх. Коригування даних 
полягає у заповненні пропусків та зменшенні викидів, що виходять за діапазон допус- 
тимих значень змінних. Нормування даних означає їх логарифмування або приведення 
до зручного діапазону їх зміни, наприклад, від 0 до І, або від - І до "1. Дуже часто для 
побудови моделі використовують не самі абсолютні значення статистичних даних, а 
їх перші чи другі різниці. Поширеним методом нормування даних є їх логарифмуван- 
ня з наступним формуванням додаткових часових рядів з перших чи других різниць. 
Як правило, із значень ряду віднімається його середнє значення, щоб отримати можли- 
вість працювати з відхиленнями, а не абсолютними значеннями змінних. Застосування 
того чи іншого методу для нормування даних може визначатися в кожному випадку 
по-своєму. 

АКФ та ЧАКФ використовують для визначення попередньої оцінки порядку мо- 
делі, тобто скільки затриманих в часі значень необхідно брати для описання процесу. 
При цьому необхідно врахувати, що АКФ дає менш «чітку» оцінку порядку процесу. 
Наприклад, для процесу АР(1) значення основної змінної у(К) та у(К - 2) будуть корельо- 
ваними, не дивлячись на те, що у(К - 2) не присутнє в моделі. Кореляція між у(К) 1 у(К-- 2), 
тобто рг, дорівнює коефіцієнту кореляції між значеннями у(К) 1 у(К - 1), помноженому 
на коефіцієнт кореляції між у(К - 1) 1 у(К - 2) або ро - рір; - р. Подібні «непрямі» ко- 
реляції присутні в АКФ будь-якого процесу авторегресії. Нагадаємо, що вибіркова АКФ 
обчислюється за виразом: 

М 
У ГО - аїїУЄК - з) - м) 


ш- Казчі 


Ре 5 


УРрюо-яї 


де У - кількість значень у вибірці даних; и - середнє значення ряду. 
Для стаціонарного процесу (це процес із постійними середнім значенням, диспер- 
сією та коваріацією) коефіцієнти г, мають нормальне розподілення та нульове середнє. 
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На відміну від АКФ, часткова АКФ між значеннями у(К) та у(К - 5) виключає вплив 
величин У(К - 1)... У(К - 5), а це означає, що коефіцієнти ЧАКФ більш чітко відобра- 
жають зв'язок між окремими значеннями основної змінної. Так, для процесу АР(1) ЧАКФ 
між у(Ю) та у(К - 2) дорівнює нулю за визначенням, що підтверджується обчисленими 
значеннями ЧАКФ. Для того, щоб знайти попередню оцінку порядку моделі, вибіркові 
коефіцієнти ЧАКФ (тобто коефіцієнти, знайдені за вибіркою даних) можна обчислити 
за допомогою простого методу, який полягає в наступному. 

а) Формують додатковий часовий ряд із відхилень основної змінної: 

О"(Юз -ДУ(Ю) -м, 
де и - середнє значення ряду. 
6) Формують рівняння першого порядку 
УЮ -Ф,, у"(К -І)че(Ю, 
де е(К) - похибка моделі. В такому рівнянні Ф!| відіграє роль коефіцієнта АКФ та 
ЧАКФ між у(Ю) та У(К - 1). Для оцінювання двох коефіцієнтів можна сформувати рів- 
няння другого порядку: 


УЮ ЕФ, У -1)-Ф,, У"(К - 2)--е(Ю), 


де Ф2г - коефіцієнт ЧАКФ між у(К) та у(к - 2). 

У загальному випадку коефіцієнти ЧАКФ стаціонарного процесу АРСС(р, 4) по- 
винні збігатися до нуля, починаючи з р-го значення. АКФ процесу АРСС(р, 4) почи- 
нає збігатися до нуля при значеннях зміщення 5 2 4. 

На другому етапі оцінюють коефіцієнти (параметри) різницевого рівняння, ви- 
користовуючи принцип економії або збереження. Цей принцип означає, що кількість 
коефіцієнтів, що оцінюються, не повинна перевищувати їх необхідне число. 

При моделюванні процесів необхідно пам'ятати, що поведінку процесу необхідно 
апроксимувати за допомогою рівнянь, а не намагатися описати її до найменших дріб- 
ниць. Необхідно враховувати також, що різні за формою моделі можуть мати однако- 
ві властивості. 

На третьому етапі отримана модель діагностується, тобто виконується перевір- 
ка на адекватність. Діагностика складається з наступних кроків. 

а) Візуальне дослідження графіка похибок моделі е(к) - у"(К) - у(К), де у(Ю - 
оцінка змінної, отримана за допомогою рівняння. На графіку не повинно бути вики- 
дів та довгих інтервалів, на яких похибка приймає великі значення (тобто довгих ін- 
тервалів значної неадекватності). 

6) Похибки моделі не повинні бути корельовані між собою. Для аналізу наявності 
кореляції між значеннями похибок необхідно обчислити АФ та ЧАКФ для ряду /е(К)/ і 
за допомогою О-статистики визначити ступінь корельованості (наприклад, О-статис- 
тика вважається несуттєвою до рівня 109). 

в) Для моделі 2 - 3 порядку оцінки параметрів повинні збігатися до усталених 
значень після 40 -- 60 ітерацій алгоритму оцінювання. Якщо кількість ітерацій набагато 
перевищує вказані числа, то це свідчить про те, що процес може бути нестаціонарним. 

г) Сума квадратів похибок повинна бути мінімальною у порівнянні з усіма інши- 
ми моделями, тобто: 


Уе(еУТУЧ(Ю- ЮР -зтіп- 


д) Для оцінки адекватності моделі також використовують інформаційний крите- 
рій Акайке 


ІКА- мі У, е «одн 
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та критерій Байєса-Шварца 
М 
КВУ мів ЗРО нти) ; 
Ка 


де п - рч411 - число параметрів моделі, які оцінюються за допомогою статистичних 
даних (р - число параметрів авторегресійної частини моделі; 4 - число параметрів 
ковзного середнього; І з'являється тоді, коли оцінюється зміщення, тобто ад). 

У правій частині виразів знаходиться сума квадратів похибок, а тому за цими 
критеріями вибирають ту модель, для якої критерії приймають найменші значення. 
Введення нового регресора приводить до збільшення критерію (при цьому збільшує- 
ться п), але разом з тим зменшується сума квадратів похибок 1 критерій в цілому змен- 
шується. Якщо регресор не покращує модель, то критерій збільшується. Необхідно 
також відзначити, що асимптотичні властивості для довгих вибірок кращі у критерії 
Байєса-Шварца. 

е) Окрім згаданих параметрів, для визначення адекватності моделі використовують 
коефіцієнт множинної детермінації А", Е-статистику Фішера та статистику Дарбіна- 
Уотсона (для перевірки корельованості похибок). 

Коректне використання методики Бокса-Дженкінса забезпечує побудову адек- 
ватної математичної моделі процесу, якщо експериментальні дані відповідають ви- 
могам представництва та інформативності. Перша вимога означає, що вибірка даних 
повинна охоплювати досить довгий проміжок часу, щоб перекривати ті режими функ- 
ціонування процесу, які цікаві для дослідника. Вимога інформативності означає, що 
вибірка повинна вміщувати в собі кількість інформації, достатню для оцінювання ко- 
ефіцієнтів моделі. Наприклад, якщо моделюється процес другого порядку, то вибірка 
повинна забезпечувати обчислення першої та другої похідної. Умову інформативності 
ще називають умовою достатнього збудження процесу. 


Ідентифікація моделей АКІМА-АРПКС зсувних 
процесів Південного берега Криму 


Ключове завдання аналізу зсувних процесів Південного берега Криму полягає в 
попередній оцінці і подальшому якісному прогнозі зсувної активності та об'єму кош- 
тів для попередження або подолання катастрофічних наслідків цих процесів. Заздалегідь 
сплановане укріплення зсувонебезпечних ділянок дозволяє раціонально розподілити 
ресурси, які виділяються на протизсувні роботи в регіоні, не зачіпаючи при цьому ін- 
тересів населення та відпочивальників на Південному березі Криму. Подібна політика 
приводить до виключення катастрофічних ситуацій, які загрожують життю та госпо- 
дарській діяльності, а також до підвищення рейтингу регіону, який використовується 
з рекреаційною, оздоровчою та туристичною метою, завдяки повній відповідності по- 
требам споживачів цих послуг. 

Організація даних 

На Південному березі Криму відповідні структури ведуть спостереження за зсув- 
ними процесами та накопичують дані, які їх описують. Загальноприйнято досліджу- 
вати відсоток зсувів, що активізувалися впродовж року, проте в моделях, що будуть 
побудовані, є сенс враховувати не відсоток, а загальну кількість активних зсувів. Це по- 
в'язано з тим, що кожного року появляються нові зсуви і вносяться до кадастру, а ста- 
рі (чи активні, чи стабілізовані) також впливають на загальну кількість, тому відносна 
кількість зсувів, врахована у відсотках, не зовсім об'єктивно описує ці процеси. Роз- 
глянемо дані з кількості активних зсувів на Південному березі Криму за кожен рік за 
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період з 1962 по 2005 роки. Дані були представлені в таблиці Ехсеї і експортовані в 
ЗТАТІЗТІСА. Таким чином, необхідно провести аналіз однієї змінної з 44 спостере- 
женнями. 

Попередній аналіз 

Спочатку побудуємо початковий ряд графічним методом (|З, а також гістограму 
розподілу накопичених даних (рис. 1). 


Гистограмма (Зриеацізічені? 10у"44с) 


Линейное вичерчиваниє (бргеацізеєіо 10м"44с) 
Оползни з 43"20"погтаїх, 134,7442, 87.827) 


Оползни 
мо огор5 


год Оползни 


У Линеднов вьнермивание Греафачені? Тис) Ї ВУ Гистогравача ГУриеакіінесі Т00с Ї 


Рисунок І - Лінійний графік і гістограма кількості активних зсувів 


З графіка змінної чітко бачимо, що існує сезонна залежність кількості активних 
зсувів. Гістограма ілюструє той факт, що розподіл ряду не є нормальним, але близь- 
кий до нього, і має сенс спробувати встановити залежності для цього ряду. 

Аналіз даних 


Знаходження залежності в представлених даних являє собою завдання: розбити 
початковий ряд на 2 складові - детерміновану функцію і чисто випадкову складову. 
Випадкова складова повинна бути рядом Гауса з незалежними прирощеннями. 

Для встановлення характеру невипадкової складової побудуємо автокореляцій- 
ну та часткову автокореляційну функції початкових даних (рис. 2). 


а-гдівв вапепагев) 


щ 101 1625 | 

15 093 1625 З 

Граница доверит. 0 Граница дорерит 
10 05 00 05 2 


ІД Резмілмлосогів Гитоїюо пола "дане До Ропаімлосоюейіог Гог | «1 
Рисунок 2 - Автокореляційна функція (АКФ) і часткова 
автокореляційна функція (ЧАКФ) 

Автокореляційна функція має тенденцію до згасання, що говорить про стаціонар- 
ність процесу. Стаціонарні ряди мають постійні за часом середнє, дисперсію й автоко- 
реляції (тобто сезонні залежності видаляються за допомогою різниць |41). Для моделі 
АРПКС необхідно, щоб ряд був стаціонарним, це означає, що його середнє постійно, 
а вибіркова дисперсія й автокореляція не міняються з часом. 

За виглядом автокореляційної функції 1 часткової автокореляційної функції мож- 
на припустити, що ряд описується моделлю авторегресії 1-го порядку. 
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Розглянемо декілька моделей аналізу часових рядів АРПКС (авторегресії та про- 
інтегрованого ковзного середнього): 

- авторегресія 1-го порядку АР(1) або АРПКС(1) 

Утефо 9 ФІЧУрб 

де ЇУ;, У; - кількість активних зсувів у відповідному році; 

- авторегресія 4-го порядку АР(4) або АРПКСІ(4) 

Уре фу РФТКореФІ Тр ФУТоо ОО 
- авторегресія 1-го порядку і ковзного середнього 1-го порядку АРПКС(І, 1) 
У; - фо РІ Тр рої 

де а; 1- білий шум; 

- авторегресія 2-го порядку і ковзного середнього 1-го порядку АРПКСІ(2,1) 

Трефо Я фІТУ; р Ж фо" Тр 2 З бз"аі- Б 
- авторегресія 1-го порядку 1 ковзного середнього 2-го порядку АРПКС(1,2) 
Уг-бо КЕ БІТ; 1 Ж БзЧаї 1 Я БзЗаї-2. 

Аналіз якості моделей проведемо за такими критеріями|3 1: 

1) коефіцієнт детермінації; 

2) середній квадрат залишків; 

3) статистика Дарбіна- Уотсона. 

Аналіз якості прогнозу проведемо за такими критеріями: 

1) середнє квадратичне відхилення - СКВ: 


1 - 
ВМУЕ СО» - у) 
ізі 


2) середня абсолютна похибка - САВ: 
І ли 
МАЕ оі3, р і 

і-і 


де у: - значення періоду, у; - прогнозна величина для і - 1,2,..., п; 


2 


3) середня абсолютна відсоткова похибка - САВП: 


МАРЕЗ ТУ РОМ 
У 


«1009; 


п Ш Дуо 
пи -У;) 


іх 2 ьо 
ФА є (У і ) Ра У (У і ) 
Поіа Па 
Результати оцінювання якості моделі та якості прогнозу для запропонованих мо- 


делей АРПКС наведені в табл. 1. 
Таблиця 1 - Аналіз якості моделі та якості прогнозу 


П/Р 


4) коефіцієнт нерівності Тейла Ю: У з 


3 Якість моделі Якість прогноз 
Мода АНИКО аа Рано З бу | скв | САВ АВ . К-т Тейла 0 
АРПКСИ) 0,501 1625 2,020 | 4032 31,39 28,895 0,139 
АРПКС(4) 0,555 1557 1,985 | 40,04 32,82 29.49, 0,137 
АРПКСИ,1) 0,512 1662 1,983 | 40,75 31,86 29,39 0,140 
АРПКС(2,1) 0,523 1625 2,022 | 4031 31,2 29,29 0,138 
АРПКСИ,2) 0,546 1583 2,053 | 39,79 32,07 28,59 0,136 
АРПКС ЩІ,2) 0,595 1528 1,960 | 48,71 32,01 26,09 0,135 


З наведених оцінок неможливо обрати модель, яка би значно відрізнялась від ін- 
ших за однією чи декількома характеристиками. Значення, наведені в стовпчиках, май- 
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же співпадають, або відрізняються зовсім мало, що говорить про однаковість моделей. 
Очевидно, що всі моделі майже еквівалентні, тобто додавання складності в модель не 
привело до значного поліпшення оцінок моделі або оцінок прогнозу. Таким чином, можна 
зробити висновок: з наведених моделей найкраща та модель, яка є більш простою, тобто 
авторегресія 1-го порядку - АРПКС(1). 

Всі розглянуті моделі для даного процесу назвемо умовно адекватними тому, що 
статистики, які їх характеризують, ледве досягають прийнятного значення, допусти- 
мого для подальшого прогнозування. Отже, отриману модель авторегресії 1-го порядку 
слід використовувати для прогнозування, але краще в сукупності з іншими моделями, 
які побудовано методами, що відрізняються від аналізу часових рядів АРПКС. Прогнози 
за різними моделями можна узагальнити 1 використовувати як прогноз, побудований 
за однією найбільш адекватною моделлю. Наведемо на рис. 3 результати прогнозуван- 
ня за допомогою моделі АРПКС(І): прогнозне значення (Еіїйеай), величину, що спос- 
терігалася, (Асішаї) та залишки (Везідшаї), а також автокореляційну функцію (АКФ) і 
часткову авторегресійну функцію (ЧАКФ) залишків даної моделі. 


ни Едиаїіоп: ОМТІПЕ Міогікбїїе: ПМ(АВІМАУ" ех) 


Мемі| Ргосз | Обіесіз| Римі | Мате! Невге| Езітаїе| Рогесазі | біаіз| Везідз| 
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Рисунок 3 - Побудова авторегресійної моделі та АКФ і ЧАКФ 


Пошук інших авторегресійних моделей 


Перетворимо дані таким чином, щоб зменшити амплітуду коливання і видалити 
уклін тренда. Попереднім перетворенням логарифмування «погасимо» слабо вираже- 
ну мультиплікативну природу початкового ряду, після цього віднімемо тренд. Ці дії 
автоматизовані в системі ЗТАТІЗТІСА |З, 1 доступні в меню інші перетворення 1 гра- 
фіки (рис. 4). 

Візуальний аналіз результатів показує, що коливання випадкової величини від- 
бувається навколо «нуля» та не перевищує одиниці за амплітудою, що свідчить про ста- 
ціонарність процесу. 

Підберемо параметри моделі, вибравши модель АРПКС(2,3), тобто р - 2 - авто- 
регресія 2-го порядку 14 - 3 - ковзне середнє 3-го порядку: 

Уагіаріе: ОПОЛЗНИ : іп(х); х-4,477-,015х- 


Тгапвеіогтаїіопз: 
Модеї: (2,0,3) 
Мо. ої орв5.: 43 Іпісіаї 55- 9,0324 Кіпа1 55- 2,9724 (32,915)  М5- ,07822 
Рагатеїкег5 (р/Р5-Айгогедге55іуе, 4/05-Момуіпд ауег.); рідпіїдьс: р«.05 


рІ1 р (2) 4(1) 4 (2) 4(3) 
Оценка: 1,3994 -,5044 ,60215 -, 4332 ,83025 
5ка.Егг.: ,16513 ,16454 ,08510 ,13395 ,08008 
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Вьмерої уагіабів: ОПОЛЗНИ 
іп(х); ж-4,477- 01571 


Ра 
ЕГ 
Га) 
Е 
о 
ЕЕ 
о 


20 25 
Номера случая 


Рисунок 4 - Візуальний аналіз результатів 
Спроможні ни параметрів підсвічуються сірим кольором. 


Гр) | | | 13993881 0,165130| 8.447444 0 000000) 1.065099 

|р(2) | | | -0,5044431 0,.1564535| -3, 06590! 0.003983! (0.837533! 0 171385| 
201) | | 0502154.0085097| 7.076120 000000) 0429885) 0,774423| 
бони 40433193! 0,133945! -3 23408 0.002627! -0,704353| -0.152033| 
0830250 0080078 10.36798 0 000000 0568140 0992360 | 


резражтав вена анна заннннняо 


Рисунок 5 - Оцінки параметрів 

На рис. 5 наведено оцінки параметрів - 1 стовпчик, стандартні відхилення (2 ст.), 
є - статистика Стьюдента (3 ст.), ймовірності відхилення цих параметрів (4 стовпчик 
«р», в ньому всі значення менші за 0,005, тобто 0,59) та нижні Й верхні границі оці- 
нок при рівні довіри 9597 (5 ст.). Середній квадрат залишків становить 0,07822. 

Для перевірки адекватності моделі скористаємося візуальними методами, представ- 
леними в системі ЗТАТІЗТІСА. Як вже було сказано, випадкова складова - залишки -- 
повинні бути нормально розподілені (рис. 6). 
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Рисунок 6 - Гістограма залишків моделі АРПКСІ?; 0; 3) 
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Внаслідок того, що прирощення випадкової складової повинні бути незалежні, 
автокореляції 1 часткові автокореляції залишків не повинні виходити за допустимі ін- 


тервали (рис. 7). 
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Рисунок 7 - Автокореляційна і часткова автокореляційна функції 


для залишків АРПКС моделі 


Як видно, автокореляції 1 часткові автокореляції (рис. 7) цілком лежать в допусти- 
мих інтервалах. (У правильно підібраній моделі залишки будуть схожі на білий шум: 
у них не буде періодичних коливань, систематичного зсуву, між ними не буде сильних 


кореляцій). 


Додатковим методом аналізу адекватності моделі є графік розташування залиш- 
ків на нормальному імовірнісному папері (рис. 8). 
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Рисунок 8 - Нормальний графік для залишків моделі АРПКС (2; 0; 3) 


Точки на графіку лежать близько до прямої, що добре характеризує модель. 


Отже, модель достатньо достовірна. Подивимося, наскільки точно вона перед- 
бачає дійсні спостереження. Для цього побудуємо прогноз за першими 7590 спосте- 
режень на тих, що залишилися 2590. 75Уо спостережень складають 33 роки, 2590 - 


11 років. 


Побудувавши графік прогнозу з 34 по 44 рік, бачимо, що прогноз добре передба- 
чає дійсні спостереження, особливо в першому періоді прогнозу. За прогнозом з 44 ро- 
ку на 2590 вперед інтуїтивно видно, що модель достатньо об'єктивна, тим більше, що 
за даними останніх п'яти років спостерігалося значне зниження активізації зсувних 


процесів. 
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тет Измеренньйй ее ПРЕДСКАЗ - - а 90.00003 Бі Ротесагіз: Модей (20.31 Сезонмов заичао 2. ВУ опесагіз- Модеі 2 0131 Сезеною змимоо 2 ЕТ | 
Рисунок 9 - Графік прогнозу на 10 років назад і на 10 років вперед 
за допомогою моделі АРПКСІ2, 3) 
Порівняння і аналіз різних моделей 


Введемо позначення АРПКС Л - авторегресія з проінтегрованим ковзним серед- 
нім логарифмована. Побудуємо 1 оцінимо моделі АРПКС Л різних порядків. Резуль- 
тати оцінювання запишемо в узагальнюючу таблицю (табл. 2). 

Таблиця 2 - Моделі АРПКС прологарифмовані із видаленим трендом 


Якість моделі Якість прогнозу 
2 Уед/п р ВМЗ5Е МАЕ МАРЕ 19 
АРПКС Л(1,1) 0,547 0,095 2,030 0,309 0,244 146,09 0,387 
АРПКС Л(1,2) 0,595 0,086 2,017 0,361 0,283 13596 0,568 
АРПКС Л(1,3) 0,590 0,088 1,997 0,354 0,282 14196 0,550 
АРПКС Л(1,4) 0,586 0,089 1,960 0,359 0,285 14390 0,554 
АРПКС Л(2,3) 0,598 0,087 2,012 0,311 0,245 13496 0,487 


Спроби провести аналіз з іншими параметрами не дали такої високої адекватнос- 
ті для моделі. Прогнозування залишків для підвищення точності аналізу виявилося до- 
сить складним завданням, оскільки залишки розподілені практично рівномірно та їх 
автокореляції близькі до нуля. Виділити сезонну складову у ряді залишків не вдалося. 

Такий прогноз можна вважати середньотерміновим тому, що термін прогнозу- 
вання становить досить тривалий відрізок часу - від року до 10. Це пов'язано зі специ- 
фікою спостережень за зсувними процесами (у зв'язку з обмеженим фінансуванням 
дані збираються раз на рік). Отже, модель АРПКС дає короткотермінові прогнози, тоб- 
то дає розрахунки на один крок вперед, але прогнози можна будувати і на триваліші 
терміни при наявності попередніх даних. 

Проте впродовж року можуть відбуватися значні зміни факторів, що впливають 
на зсувні процеси Південного берега Криму, тому є сенс в побудові нової моделі, яка 
може вчасно реагувати на різні збурення факторів, тобто модель надаватиме коротко- 
терміновий прогноз, який включатиме строк від одного тижня до двох - трьох місяців. 
Цим вимогам цілком відповідає мережа довіри Байєса (51. 


Висновки 


У статті розглянуто аналіз часових рядів за допомогою АРПКС моделей зсувних 
процесів на Південному березі Криму. Наведена методика Бокса-Дженкінса для побу- 
дови математичної моделі процесу у вигляді різницевого рівняння. Розглянуто деталь- 
но етапи побудови моделі. 

Проведено ідентифікацію моделей АКІМА-АРПКС зсувних процесів Південно- 
го берега Криму 1 визначено, що при ускладненні моделі додатковими параметрами 
оцінки адекватності моделі й адекватності прогнозу не покращуються. 
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Розрахунки оцінок моделей АРПКС(1), АРПКС(ТІ,1), АРПКС(1,2), АРПКСИ,3), 
АРПКСІ(1,4) підтверджують результати попереднього аналізу. З таблиці видно, що мо- 
делі мають майже однакові статистичні характеристики у порівнянні одна з одною. 
Оцінки коефіцієнтів моделі суттєво відрізняються від нуля на рівні 195, а корені ха- 
рактеристичного рівняння знаходяться всередині кола одиничного радіуса. Значення 
О-статистики свідчать про те, що автокореляція між похибками є статистично несут- 
тєвою, тобто нуль-гіпотеза О - 0 підтверджується. Побудовано графіки залишків ав- 
токореляційної та часткової автокореляційної функцій (ЧАКФ). 

Перед оцінюванням наступної моделі проведена попередня обробка даних - ло- 
гарифмування та видалення тренду. Іншим підходом до попередньої обробки даних 
може бути використання квадратного кореня значень або іншої процедури. 

У вищенаведених моделях сезонний ефект не був описаний за допомогою лагу. 
Крім такого підходу до моделювання сезонного ефекту існують також 1 інші, які мо- 
жуть дати точніше описання процесу. Наприклад, існують моделі з мультиплікативни- 
ми сезонними компонентами 

а-адууЮ«чнвИзачВІє(Ю, 
де В8,є(К - 4) входить в модель у мультиплікативній формі. Така модель може мати кращі 


показники прогнозування. Наступним кроком дослідження процесу може бути тесту- 
вання часового ряду на наявність гетероскедастичності, тобто чи є дисперсія ряду 
змінною величиною. 
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В.Н. Таран, А.О. Пашко 

Оценка прогнозов оползневьїх процессов Южного берега Крьтма 

с помощью анализа временньгїх рядов (АРПСС) 

В статье исследованьт общиєе закономерности построения авторегрессионньх моделей оползневьіх 
процессов на примере Южного берега Крьтма. Проведена идентификация АРПСС-моделей, вьшолненьт 
расчетьт оценок параметров зтих моделей. Вьбраньт критериий оценки моделей и оценки прогнозов, по 
которьм вьіделяется модель среди остальнь(хх, а также построен прогноз на 10 лет. 
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Езбітабіоп ої Ргоспо5бе5 ої Г апдзіїде Ргосе55ез ої Ше Зошірпегп Соазі ої Сгітеа 

Бу Меапз5 ої Апаїубіз ої Тетрогаї! Коубз5 (АВІМА) 

Сепегаї сопбоглійез ут а Іаму ої соп5ігасіїоп ої аціогевтеззіме плодсеіз ої Іапд5іде ргосе55е8 аге іпус5іватед, 
оп Фе ехатріе ої Ше 5оцірегп Соа5і ої Сгітеа. Аиірепісайоп ої АКІМА-тодеія 15 сопдистіва, Ше 
саїсиіайоп5 ої ез(тайопз ої рагатеїегя ої пе5е плодеї8 аге ехесиісад. Те сгіегіа ої езітайоп ої тодеіз 
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ргоєпо5зі8 15 Бицї Гог 10 уваге5. 
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